Effects of current on vortex and transverse domain walls 
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By using the spin torque model in ferromagnets, we compare the response of vortex and transverse 
walls to the electrical current. For a defect-free sample and a small applied current, the steady state 
wall mobility is independent of the wall structure. In the presence of defects, the minimum current 
required to overcome the wall pinning potential is much smaller for the vortex wall than for the 
transverse wall. During the wall motion, the vortex wall tends to transform to the transverse wall. 
We construct a phase diagram for the wall mobility and the wall transformation driven by the 
current. 

PACS numbers: 



I. INTRODUCTION 

In a magnetic wire, there are two standard walls: 
transverse wall (TW) and vortex wall (VW). In a cer- 
tain range of the wire width and thickness, both types of 
the walls are stable [JQ. 

Experimentally, one can create 
either type of the wall by applying a magnetic field in var- 
ious directions 0, Q| . The dynamics of both walls driven 
by a magnetic field has been extensively simulated 0, HI • 
In general, the mobility of the vortex wall is smaller than 
the transverse wall, and at a large magnetic field, the 
vortex core tends to move toward the wire edge and may 
be annihilated at the sample boundary. Thus a transfor- 
mation from vortex wall to transverse wall may occur. 
In this paper, we address the wall dynamics driven by a 
current. It would be interesting to see how the current- 
driven wall dynamics qualitatively differs from the field- 
driven wall dynamics. The present work is also moti- 
vated by recent experimental observations on the domain 
wall motion £|. Particularly, Klaui el al. have quantita- 
tively measured the correlation between wall mobilities 
and wall structures, and a transformation from vortex 
wall to transverse wall has been clearly identified. 



II. MODEL 

We start with the generalized Landau-Lifshitz-Gilbert 
equation which includes the recently derived spin transfer 
torques 0: 
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where 7 is the gyromagnetic ratio, and H e // is the 
effective magnetic field including the external field, the 
anisotropy field, the magnetostatic field, and the ex- 
change field, a is the Gilbert damping parameter, bj = 
PjefJ-B I 'eM s , Cj — £bj, P is the spin polarization of the 
current, j e is the current density along the length di- 
rection of nanowire, /is is Bohr magneton, and £ is a 



dimensionless constant which describes the degree of the 
nonadiabaticity between the spin of the nonequilibrium 
conduction electrons and local magnetization. 



III. CALCULATION PROCEDURE 



We have solved the above LLG equation by performing 
a routine micromagnetic simulation for a magnetic wire 
whose geometrical size is 2/im long (x-direction, also the 
current direction), 128nm wide, and 8nm thick. The 
grid size is chosen as 4 x 4 x 8nm 3 . The magnetization 
at both ends is set to be along the x-direction and direct 
inward to the wire; and we use free boundary conditions 
for other dimensions. Since the domain wall may move 
up to several micrometers in some cases to be discussed 
below, it is important to keep the domain wall far away 
from the ends. To reduce this end effects due to the 
finite length of the wire, we shift the domain wall to the 
center of the wire after every small displacement or time 
interval, corresponding to the geometrical structure of 
wire. Using the similar initial domain configurations as 
used in Q] , these two types of walls will form and they are 
our initial wall configurations at t = 0. Starting at t = 0, 
a static spin torque or a static magnetic field is applied 
to the wall until the end of the simulation. The strength 
of the spin torque or the field is chosen in a certain range 
and the dynamics of the wall after t > is reported 
in the paper. In simulating the pinning potential, we 
choose a geometrical boundary roughness by periodically 
removing, for every 32 grids in the length direction, one 
grid from each edge in the direction of wire width 8] . 

The material parameters are: the exchange constant 
A = 1.3 x 10~ 6 erg/cm, the anisotropy field Hk — 
0(Oe), the saturation magnetization M s — SOOemu/cc, 
the damping parameter a is 0.02, the spin polarization 
P = 0.5, and the non-adiabaticity factor £ = 0.04. 
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inside the wire only when the current /field is sufficiently 
small. In this case, the wall motion remains steady and 
the velocity is given by Eq. (2) as shown in Fig. (I). The 
dynamics of the vortex wall at a higher current/field will 
be discussed in next section. For the transverse wall, if 
one applies a larger current/field, the uniform motion be- 
comes unstable. In fact, the wall will constantly deform 
during its motion and the velocity is highly non-uniform, 
known as Walker's breakdown |12j . In the film we con- 
sider here, we choose the applied field and current far 
below the critical values for Walker's breakdown. We do 
not intend to describe the complicated dynamics after 
the breakdown in this paper. 



B. Domain wall motion in the presence of defects 



FIG. 1: Domain-wall steady state velocity v x as a function 
of (a) current J e , and (b) field H e . In (a) v x = —cj/a and in 
(b) v x = W /a, both agree with the analytical result in 
(2) 



IV. RESULTS 
A. Defect-free domain wall motion 

A simplest case is a defect-free wire. In this case, the 
steady state wall velocity v x can be analytically derived, 

lH e W CJ 

v x = (2) 

a a 

where H e it the external field and W is the domain wall 
width defined as 

where S is the cross section area of the wire. The as- 
sumption used in deriving Eq. (2) is that the domain 
wall moves uniformly, i.e., M = M(r — v x tx), see pa- 
per of this proceedings. Thus, in the absence of the 
field, the steady velocity is —cj/a, independent of the 
wall structure In the presence of the field, however, 
the velocity depends on the wall width defined in Eq. (3). 
The vortex wall has a smaller effective width compared 
to the transverse wall and thus the vortex wall has lower 
velocity In Fig- CO) we show the steady state veloc- 
ity for the vortex wall and the transverse wall. The wall 
velocity is well fitted by Eq. (2). 

However, the uniform motion is expected to break 
down even for the defect- free wire when the external field 
or current density is large. For the vortex wall, when a 
current/field is turned on, the wall moves along the wire, 
and at the same time the vortex core also tends to move 
perpendicular to the wire. A restoring force [T(i| may stop 
this perpendicular movement and keep the vortex core 



While there are many types of defects to pin a domain 
wall, we restrict our discussion on the defect induced by 
the roughness on the edges of the wire width. Specifically, 
we periodically remove one grid at each edge as described 
above in otherwise a defect-free wire. 

Let us start with the transverse wall. When a current 
is applied, the domain wall begins to move and in the 
meanwhile the wall width shrinks. The wall shape re- 
mains similar to the original one. If the current is small, 
the wall eventually stops. The displacement of the center 
of the wall is proportional to the current density and the 
magnitude could be as large as several tens of nanome- 
ters. The terminal average velocity is zero. When the 
current exceeds a critical value J tl , as shown in Fig. (2a), 
the domain wall is able to overcome the pinning and can 
sustain a significant velocity, which is oscillating corre- 
sponding the defect positions (not shown here). In fact, 
the average velocity are smaller but still close to —cj/a 
once the current density is larger than the critical cur- 
rent. It is understood that part of the effect of the defects 
has been compensated by the deformation of the domain 
wall and the average velocity has not been reduced sig- 
nificantly. 

One interesting feature of the transverse wall is that 
the wall velocity as a function of the current is hysteretic. 
If the wall reaches the critical current, J tl , the wall begins 
to move with an average velocity, shown in Fig. (2a). 
Then after the wall is virtually depinned (J > J tl ) and a 
uniform wall velocity is established, one begins to reduce 
the current. It is found that the domain wall remains 
to have a finite velocity even if the current density has 
reduced to a value lower than Jt ± . One needs to further 
reduce the current density until reaching a second critical 
current, J t2 , in order to completely stop the domain wall, 
see Fig. (2a). The hysteretic feature may be qualitatively 
understood in terms of the competition between the wall 
kinetic energy and defect potential 9] . When the current 
density is reduced from a higher value, the velocity of the 
wall remains large and it would be difficult for defects to 
pin the wall. On the other hand, when we increase the 
current density from a smaller value, the wall is originally 
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FIG. 2: The dependence of average domain- wall velocity on 
the current and wall structures. In (a), solid line shows the 
velocity of transverse wall by applying a current from the 
static state (v x (0) = 0,J e (0) = 0); dashed line represents 
the velocity by applying a reduced current from a depinned 
state (v x (0) 7^ 0, J e (0) = Jti)- In (b) the open circle line 
shows the average velocity of vortex wall, which depins at J V1 
and transforms to a transverse wall at J V2 ; the up-triangle 
line represents the velocity of the transverse wall after the 
transformation: upon J V3 , the transformed wall keeps the 
same velocity as in (a) with corresponding currents. 



This is somewhat understandable because the effects of 
edge roughness on a vortex wall is quite small. In the 
previous work [l^, it has been shown that the similar 
defect may trap the domain wall when being present at 
the center of vortex core, and it has almost no effects to 
the wall if staying at the edge. 

When the current continues to increase, the average ve- 
locity of the wall proportionally increases. In the mean- 
while, the center of the vortex core moves closer to the 
edge. At a second critical J V2 , the vortex reaches the 
boundary and the vortex wall is transformed into a trans- 
verse wall. Because the critical current for the transverse 
wall Jt x is much larger, the wall motion stops. These 
features seem to agree with the experimental observa- 
tion [4j. If the applied current is further increased upon 
to the third critical value J V3 which is smaller than J tl , 
the wall can still keep a finite velocity after the trans- 
formation is completed. And this average velocity just 
equals to the terminal average velocity that the trans- 
verse wall may sustain after the current is reduced from 
higher level, i.e.Jt 1 in Fig. (2a). It is also explained as 
the hysteretic feature of the transverse wall we described 
above: the kinetic energy gained before the wall trans- 
formation helps the transformed wall sustain a significant 
velocity. 



V. CONCLUSIONS 



at zero velocity and thus the pinning is more effective. 

The motion of the vortex wall is far more complex than 
that of transverse wall, because the transformation hap- 
pens when the applied current/field is high. As in the 
transverse domain wall, one needs a critical value of the 
current to maintain a finite terminal velocity with vortex 
core inside the wall. However, the most interesting fea- 
ture is that the critical current for the vortex wall, J Vl , 
is far smaller than J tl , i.e., the vortex wall is less sensi- 
tive to the defect potential. As shown in Fig. (2b), the 
vortex wall begins to move at a current density an order 
of magnitude smaller than that for the transverse wall. 



We have investigated the domain wall dynamics for 
vortex and transverse walls by using the spin torque 
model. Interestingly, while the terminal velocity is in- 
dependent of the wall structure for a steady wall motion 
in defect-free wire, the depinning current of the vortex 
wall is an order of magnitude smaller than the transverse 
wall in the presence of defects. However, the vortex wall 
is not as stable as the transverse wall during its motion. 
The vortex core may be annihilated at the edges and the 
transformation to the transverse wall occurs. These con- 
clusions seem to explain the experimental observation Q . 
This work is supported by NSF-DMR-0314456 and the 
MU research board. 
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